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Recent developments in analog computers have aroused fresh in- 
terest in automatic iterative computation. While the technique is 
not new, it has not been widely used with analog computers in the 
past, and is therefore unfamiliar to many analog programmers. It 
is the purpose of this report to provide an introduction to the tech- 
niques involved and to suggest representative applications. 

While the analog computer is, in many ways, well- suited for this 
mode of operation, many of the basic concepts (logical decision, 
storage, discrete variables) are more familiar to users of digital 
computers. Iterative analog computation is, in fact, a form of hy- 
brid computation, employing both analog and digital techniques, 
and large and complex iterative problems are most effectively at- 
tacked by combining digital devices (flip-flops, shift registers, 
gates, etc.) with conventional analog components. Nevertheless, a 
good many problems can be solved with analog components alone, 
and this report will deal almost exclusively with this type of prob- 
lem. Specifically, the report will be concerned with what can be 
done with existing components on a medium sized general purpose 
analog computer such as the PACE® TR-48. 



Of particular interest in this regard, however, is the hybrid imple- 
mentation (described in Appendix 3) of the illustrative analog 
problem solved in the body of the study. Using the EAI DES-30 
Digital Expansion System ... a general purpose, patchable digital 
logic package ... to provide digital expansion of the TR-48 analog 
computer, the all-analog simulation of this one-parameter optimizer 
or boundary value problem is modified so that the necessary logic 
functions are performed by the synchronous logic elements of this 
completely self-contained unit. 

Although capable of operating autonomously as an aid to digital 
instruction or design, the DES-30 digital logic system provides 
basic hybrid capabilities to the small computer facility. It is 
readily combined also with other general purpose analog or digital 
computers to provide extended control functions. 
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AUTOMATIC ITERATIVE OPERATION 
ON AN ANALOG COMPUTER 



by 

George Hannauer, 
Omri Serlin, and 
Peter J. Holsberg 



1. DEFINITION OF ITERATIVE OPERATION 



By iterative operation of an analog computer will be 
meant any sequence of operations that involves: 
(1) automatic cycling of the computer mode so that 
the computer goes through a number of OPERATE 
cycles, and (2) automatically changing some con- 
dition of the problem — an input, a parameter 
value, or an initial condition — between successive 
OPERATE cycles so that a slightly different prob- 
lem is solved with each iterationi'? • 

The second condition is important. It distinguishes 
true iterative operation from repetitive operation, 
which satisfies only the first condition. Repetitive 
operation, in fact, is not a computational technique 
at all, but merely a form of readout, or display. No 
special techniques are involved in programming a 
rep -op problem; in fact, programming a problem 
is the same for rep-op as for ''real time" read- 
out, with the exception that the programmer must 
restrict himself to high-speed components (e.g. 
quarter- square multipliers instead of servos). 

In repetitive operation, much of the information 
generated is redundant; except for occasional 



manual parameter changes, the same problem is 
solved again and again, and the same curve gen- 
erated repeatedly. In contrast, each OPERATE 
cycle in iterative computation investigates a slightly 
different set of parameter values, and produces 
slightly different results. In a given number of 
cycles, therefore, iterative operation generates 
considerably more information, utilizing more fully 
and efficiently the analog computer's ability to 
produce a large amount of data in a short period 
of time. 



The word ''automatic" is also an important dis- 
tinction in the above definition. In a loose sense, 
almost all operations of an analog computer may 
be classed as iterative, since rarely if ever is an 
operator satisfied with a single run. More typically, 
he will make many runs, varying parameter values 
from one run to the next, and observing the effect 
on the solution. What distinguishes iterative com- 
putation, as defined above, is the possibility of 
automating this procedure. 
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2. TYPICAL APPLICATIONS 



a) Parameter Sweep . For a simple example of 
automatic iterative operation, consider the ef- 
fect of an adjustable parameter, « ^ on an analog 
computer solution. It might be desirable to 
study the steady- state error of a position 
servo as a function of controller gain or some 
other control parameter, or the area of in- 
terest might be in the miss distance of a mis- 
sile as a function of initial launching error. 
In either case, prime interest is in a single 
number, X, (steady-state error, miss distance, 
etc) which is computed during an OPERATE 
cycle using a particular value of some param- 
eter, « (controller gain, launching error, etc.). 
A value of '« is chosen at the beginning of the 
OPERATE cycle, and the corresponding value 
of X is computed (usually X will be a final 
value, i.e. the output of some component at the 
end of the OPERATE cycle). It is desired to 
obtain a plot of X versus « . 

Obtaining such a plot manually can be quite 
time-consuming. The operator has to choose 
several values of « , changing pots manually 
between runs and tabulating the result X after 
each run. He can then make a point-plot of X 
versus « and connect the points by a smooth 
curve. The procedure outlined below replaces 
this tedious work with a single computer run, 
in which the desired graph is automatically 
produced by an X- Y plotter. The necessary 
circuit is described in Appendix 1. 

b) Optimization . Optimization problems arise 
quite frequently in design studies. In many 
cases, a long sequence of computer runs is 
made, with parameter adjustments between 
runs to achieve some optimum condition such 
as maximum profit, minimum cost, minimiun 
error, etc. 

For example, consider the design of a control 
system to stabilize the orientation of a satellite. 
Suppose the system is characterized by two 
parameters: a gain IC and a dead- zone d,. The 
problem is to correct a given orientation error 
with minimum fuel consumption. The general 
approach is given in Figure 2-1. 

The computer is placed in the OPERATE mode 
for a sufficient length of time to correct the 
initial error, and the amount of fuel used is 
remembered" by a storage device. The values 
of X from several successive OPERATE cycles 
may be stored if needed. On the basis of these 



♦ CIRCUIT FOR SOLVING 
EQUATIONS OF MOTION 
d ^ OF THE SATELLITE 



CONTROL CIRCUITRY 


1) 


RESET COMPUTER, BUT STORE 




THE VALUE OF X OBTAINED IN 




COMPUTATION. 


2) 


CALCULATE CHANGES IN K 




AND/OR d TO DECREASE X. 


3) 


OPERATE COMPUTER WITH 




NEW PARAMETER VALUES. 


4) 


REPEAT STEPS 1-3. 



d = WIDTH OF DEADZONE IN ERROR SENSING DEVICE 
K = CONTROLLER GAIN 

X = AMOUNT OF FUEL USED IN THE CORRECTING 
MANEUVER 

Fiaure 2-1 

results, changes are made in d or K or both, 
and the problem is re-run. The process is 
continued until a minimum is reached. 

This type of problem differs from parameter 
sweep in that the former is purely '*open- 
loop" i.e. the changes in parameter values 
are determined in advance, while in optimiza- 
tion problems the parameter changes are 
allowed to depend on the results of previous 
runs. The circuitry for performing the re- 
quired storage, parameter changes, etc. will 
be discussed in Sections 3 and 4, 

c) Boundary- Value Problems . The analog com- 
puter is essentially an initial-value device; that 
is, from the many possible solutions to a given 
set of differential equations, a single one is de- 
termined by specifying initial values for the 
variables, i.e. values at time zero. In contrast, 
many problems are stated in terms of boundary 
values, which must be met by trial and error, 
i.e. by some sort of iteration nrocess. 

The most common examples of boundary value 
problems arise in the solution of partial dif- 
ferential equations which are discussed below. 
However, occasionally one runs into the same 
problem with ordinary differential equations, 
as in the following example: 

A bullet is fired from a cannon at an inclina- 
tion e from the horizontal, and with a muzzle 
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velocity Vq, In the absence of air resistance, suggest that the decay should satisfy the dif- 

the trajectory of the bullet is described by the ferential equation. 

equations: 

C=-kCn (2-3) 

X =Vo cose x(o) =0 (2-1) 



y = -g y(o) = 0; y(o) = Vq sin 6 (2-2) 

where x and y are the horizontal and vertical 
coordinates of the bullet respectively, as in 
Figure 2-2. 



y 




Figure 2-2 

The problem is to find the firing angle Q that 
will make the bullet strike a target at ground 
level, at a distance R from the point of firing. 
Thus, the problem is to find the value of x(t) 
when y(t) becomes zero. Figure 2-3 gives a 
general flow diagram for solving this type of 
problem. Sections 3 and 4 contain descriptions 
of the necessary circuits. Section 5 contains 
the solution of this problem. 





ANALOG CIRCUIT FOR SOLVING 
EQUATIONS 2-1 AND 2-2 


X(T) 












CONTROL CIRCUITRY 

1) STORE THE VALUE OF X(T) WHEN 
Y(T) CROSSES ZERO. THIS IS 
THE RANGE R 

2) RESET COMPUTER AND CALCULATE 
A NEW VALUE OF 6 TO BRING R 
CLOSER TO THE DESIRED VALUE. 

3) OPERATE COMPUTER WITH THIS 
NEW VALUE OF 0. 

4) REPEAT STEPS 1-3. 


X(T) 







Figure 2-3 



d) Curve- Fitting. In many applications, a math- 
ematical model with free parameters must be 
fitted to empirical data, i.e. a theoretical 
curve must be passed through a number of 
experimental points. As an example, suppose 
bench- scale experiments provide data on the 
decomposition of a component in a chemical 
reaction. Chemical theory and/or experience 



where C (t) is the concentration of the chemi- 
cal at time t. The problem is to find numeri- 
cal values of k and n so as to provide the 
closest agreement between the theoretical 
curve and the experimental data. 

Let C(t) be the solution to the equation (with 
given values of k and n) and C*(t) be the ex- 
perimental curve. Then, it is desired to 
choose k and n so as to make C(t) and C*(t) 
as close as possible. 

To automate this, it is necessary to define 
some measure of the error. Defining 

e = y [c (t) - C* (t)] ^dt (2-4) 

the computer can be programmed to make e as 
small as possible. Thus, curve-fitting prob- 
lems can be treated as a special case of opti- 
mization problems. 

6) Serial Solution of Partial Differential Equa- 
tions. The solution of partial differential 
equations on a general -purpose analog com- 
puter by the usual finite-difference techniques 
involves a large number of similar circuits 
operating simultaneously. If the mesh size is 
small or the number of independent variables 
is large, then the duplication of equipment may 
make this approach prohibitively expensive. 
The simularity of circuits leads to the idea 
of using a single * 'time- shared'' circuit; that 
is, a single circuit which performs sequen- 
tially the calculations that would otherwise- be 
performed simultaneously by many similar 
circuits. 

Consider the heat equation 

= K — (2-5) 
, 2 Bt 
Bx 

with boundary conditions 

T (o, t) =f (t) 
T (L, t) =g(t) 
T (X, o) =h (X) 

This system describes the temperature T(x, t) 
at position x and time t in a thin, homogeneous 
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rod of length L with an initial temperature 
distribution h(x) and time- varying tempera- 
tures f(t) and g(t) established at the ends. 

The usual approach to analog solution of this 
problem is to treat x as a discrete variable, 
replace partial derivatives with respect to x 
by finite-difference approximations, and solve 
the resulting ordinary differential equations 
in the conventional manner. 

For serial solution, the best- approach is to treat 
t as a discrete variable and work with a system 
of ordinary differential equations in x. ^lf^:tHF 
-j:!eader- doubts4his; he"m 

with disGrete X affl t and struggle witli 

.the boundary conditions.) 

Replacing the time-derivatives with finite- 
difference approximations, the following equation 
is obtained: 



(X, t) + K 



T (x, t - At) 



At 



(2-6) 



computer during a single OPERATE cycle is an 
instantaneous temperature profile. Establishing 
boundary values T (o,t) and T (L,t) for T at both 
ends of the rod means establishing initial and final 
values for the output of integrator 2. 

Establishing the initial value is no problem since 
this means putting an initial condition directly on 
integrator #2. To establish a final value, an ap- 
propriate initial condition must be picked for 
integrator #1. This requires an iterative procedure 
similar to the one used for the missile trajectory 
problem described above. The sequence of steps 
would be as follows: 

1) Use the circuit in Figure 2-4 to calculate 
T(x, At) with the given function h(x) = 
T(x, o) as an input. Use f(At) = T(o, At) as 
an initial condition on amplifier 2, and 
guess at an initial condition for amplifier 1. 

2) At the end of the OPERATE cycle, com- 
pare the final value T (L, At) with the de- 
sired value g(At), and change the IC on 
amplifier #1 accordingly. 



where T , T refer to derivates with respect to 
X (not t). Since the integrators on an analog com- 
puter integrate with respect to time, it is necessary 
to represent distance in the heat equation by time 
on the machine. The circuit for solving equation 2- 6 
is given in Figure 2-4. 



-T(x,t-At) FROM STORAGE 




TO 

STORAGE 



Figure 2-4. Circuit for Solvinq Equation 2-6 

The fact that computer time represents distance in 
the original equation must be taken into account 
in interpreting the results, -'^'■■^"'■^k. v K\i ; 

The length of the OPERATE cycle is proportional 
to the length of the rod. The graph produced by the 



3) Iterate until the condition T(L, At) =g(At) 
is met. 

4) When this condition is met, store the com- 
plete curve T(x, At) and repeat steps 1-3 
with T(x, At) as an input to generate 
T(x, 2At). 

Continuing in this manner, a sequence of temp- 
erature profiles is obtained: T(x, o) (given), 
T(x, At), T(x, 2At). . . .T(x, nAt). 

NOTE that this application involves an iteration 
within an iteration. Such computations are referred 
to as nested computations. 



All of these problems, with the possible exception 
of the last, can be solved on the TR-48 provided 
the number of adjustableparameters involved is not 
too great. What distinguishes the last example is 
the need for curve storage — the entire curve 
T(x,At) must be stored and used as an input to 
generate T(x, t + At) in the next OPERATE cycle. 
In contrast, only a few numbers need be remem- 
bered at any one time to solve the other problems. 

Curve storage is possible on the TR-48 but the 
equipment requirements are severe. At present, it 
is more practical to solve problems requiring 
curve storage on a larger machine such as the 
100 amplifier 231-R. xA number of circuits are 
available, (References 2, 3, and 4). ' • . ^ 



3. ITERATIVE REQUIREMENTS AND CIRCUITS 
USING ONLY ANALOG INTEGRATORS 



In each of the preceding applications the conven- 
tional analog circuit that solves the dynamic equa- 
tions of the basic problem is supplemented by a 
number of additional components to perform the 
iteration. These additional components must be 
capable of performing the following functions: 

Storage; The computer must remember" 
the results of previous computations, using 
them, if necessary, in the computation of new 
parameter values for successive runs. 

Logic; The computer must decide which pa- 
rameter(s) to alter, and must route signals 
into and out of storage in response to the re- 
sults of past computations. 

Automatic Parameter Adjustment ; The com- 
puter must be able to change the values of 
important parameters and initial conditions. 

Mode Control : The computer must be auto- 
matically cycled through a sequence of modes. 
Normally, the RESET and OPERATE modes 
will alternate, but more complicated sequences 
are sometimes desirable. 

The ability of a computer to store information, 
make logical decisions, and alter its own program 
has long been a familiar concept to digital com- 
puter users. What is not so well known is that an 
analog computer can also perform these functions, 
as the following pages will indicate. 

This section will deal with the above requirements 
and how they can be met on an analog computer. 
The circuits given below are not necessarily the 
best or most economical. They have been chosen 
with one overriding consideration in mind — they 
can be mechanized by ordinary general -purpose 
analog components, without special -purpose ex- 
ternal equipment. 

a) Storage. Since an analog signal is a voltage, 
it may be stored on a charged capacitor. Read- 
out of the stored information must take place 
without drawing appreciable current from the 
capacitor, as drawing current will partially 
discharge it, resulting in an erroneous read- 
ing. This requirement for '^nondestructive 
readout" implies that an amplifier will be 
necessary. 

Consider the circuit in Figure 3-1 in which an 
amplifier and a capacitor are combined. Were it 



not for the presence of the capacitor, this circuit 
would simply be an inverter with the relation 

^out = "^in- 

The effect of the capacitor is to make Vout l^g 
slightly behind -Vin, but this lag will be negligible 
provided V^n is not changing too rapidly, and the 
capacitor is sufficiently small. (The transfer func- 
tion of the circuit is 

Vout _ 1 

and this is nearly -1 as long as the time constant 
RC is small.) 

If the relay is opened, the amplifier behaves like 
an integrator in the HOLD mode, and Vq^^. remains 
constant at the value it had when the relay opened. 
The amplifier is said to track , or follow its input 
voltage when the relay is closed, and store , or 
hold its value when the relay is opened. 




Figure 3-1 



The circuit in Figure 3-1 is essentially the initial 
condition circuit of an integrator. If an input vari- 
able is patched into the IC terminal of an integra- 
tor, the integrator will track that variable in the 
RESET mode and store it in the HOLD mode (and 
even in the OPERATE mode if the integrator has 
no other inputs) . 

Since switching an integrator from tracking to 
storing is accomplished by changing its mode, it 
is necessary to be able to control the integrator 
modes individually. This can be done ontheTR-48, 
as the integrator relay coils appear on the patch 
panel. A detailed description of TR-48 mode con- 
trol is given in Appendix 2. The devices described 
below all use the track- store principle outlined 
above but they differ in their mode control. 
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INPUT 

^ ^^^X ^ OUTPUT 

TRACKING AMPLIFIER 
Patch RESET coil to OPERATE bus 
Patch OPERATE coil to RESET bus 
Remove jS plug 

(Tracks when computer is in OPERATE; stores 
when computer is in RESET) 



Notes: 

1) The T and S amplifiers are named after 
what they do in OPERATE — T tracks 
and S stores. In RESET, they do the 
opposite. 

2) Input resistor networks need not be 
patched. 

3) Removal of plug reduces the tracking 
error by shortening the time constant. 



INPUT 




STORAGE AMPLIFIER 

Patch OPERATE and RESET coils normally 
Remove j3 plug 

(Stores when computer is in OPERATE; tracks 
when computer is in RESET) 



In order to store a result from a previous compu- 
tation cycle, a tracking amplifier and a storage 
amplifier in cascade are used, as in Figure 3-2. 
The problem variable, X (t), is tracked during an 
OPERATE cycle, and its final value is stored 
between OPERATE cycles.. The result is shown in 
Table 3-1; after every OPERATE cycle except the 
first, the output of amplifier S is the final value of 
X from the previous computer run. If it is desired 
to store a value of X other than a final value, the 
T amplifier can be replaced with a comparator- 
controlled track- store unit. 



INPUT 




OUTPUT 



COMPARATOR CONTROLLED TRACK-STORE 
UNIT 

Omit mode control connections 

Connect SJ^ (not SJ) to B through comparator 

contacts 

(Stores when comparator is open, tracks when 
comparator is closed) 
Remove j3 plug 



In OPERATE, T tracks and S stores. In RESET, 
their functions are reversed. 
x(t) 




Figure 3-2 




Table 3-1 indicates the behavior of T and S ampli- 
fiers, where X]^, X2, X3 are the solutions produced 
during the 1st, 2nd, 3rd OPERATE cycle. 



TABLE 3-1 



Computer 
Mode 



Behavior of 
Amplifier T 



Behavior of 
Amplifier S 



R 


Stores (output = 0) 


Tracks (output = 0) 


0 


Tracks - X^ 


Stores (output is still 0) 


R 


Stores - (final value) 


Tracks X-|^ (final value) 


0 


Tracks - X2 


Stores X^ (final value) 


R 


Stores - X2 (final value) 


Tracks X2 (final value) 


0 


Tracks - X 


Stores X^ (final value) 



b) Relay Logic. Logical decisions on the analog 
computer are often made by relay compara- 
tors. The relay comparator consists of a high- 
gain amplifier whose output drives the coil of 
a DPDT relay. The relay contacts normally 
make to the side marked they transfer 

to the side marked when the sum of the 
input voltages becomes negative. (See Figure 
3-3.) 




SIMPLIFIED SCHEMATIC 



Figure 3-3. The Relay Comparator 

It should be noted that the input resistors on the 
TR-48 comparator are not equal. Therefore, a pot 
may be necessary to attenuate one of the inputs 
and assure that the switching occurs exactly when 
X + Y crosses zero. 

The comparator is numbered with the number of the 
amplifier inmiediately to its left on the patch panel; 
this convention makes it easy to locate the com- 
parator when patching. The letters *'U'' and *«L" 
on circuit diagrams refer to the upper and lower 
sets of contacts respectively. 

c) Automatic Parameter Adjustment . A param- 
eter in an analog computer circuit is a num- 
ber which is constant during any single OP- 
ERATE cycle, but which is allowed to change 
value between successive OPERATE cycles. 
Usually such parameters are represented by 
pot- settings, and the pots are adjusted man- 
ually between computations. However, in auto- 
matic iterative computation, the parameter 
value is changed by the computer itself, which 
means that the parameter must be regarded 
as a computer variable, that is, as an ampli- 
fier output voltage. 



Let « be the parameter that must be maintained 
constant during a computer run and vary between 
runs by an amount A o.. The value of A « i.e. the 
magnitude and direction of the change in « ^ can 
be made to depend on the results of previous com- 
putations. Hence, A « will also be an amplifier out- 
put voltage. Methods of calculating A « are dis- 
cussed in Section 4. 

The easiest way to change « in the desired manner 
is to use an integrator, as in Figure 3-4. 



±10 




Figure 3-4 

Since « is to be constant during a computer run, 
the integrator must be in the HOLD mode when the 
rest of the computer is in the OPERATE mode. 
Since « is supposed to change between computer 
runs, the integrator must be in the OPERATE mode 
when the rest of the computer is in RESET and 
the problem is being reset for the next run. This 
type of mode behavior can be achieved on a TR-48 
by patching the OPERATE coil to the RESET bus. 
(See Appendix 2 for an explanation of TR-48 mode 
control.) 

If the computer resets for T seconds between 
OPERATE cycles, the integrator in Figure 3-4 
will integrate the constant -A^/T for T seconds, 
and its output will change by an amount Aocfor the 
next computation. 

An initial condition" must be chosen for oc, that 
is a value for the first run. This is, of course, an 
initial condition on the integrate^. The mode control 
must be arranged so that the integrator is in the 
RESET mode at the beginning of the iteration. A 
method of accomplishing this is given in Section 3d. 

The three most common uses of the signal are: 

1) as an initial condition 

2) as an additive constant, or bias voltage 

3) as a coefficient 
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In cases 1) and 2) the voltage « may be fed directly 
to the problem, as a voltage input is what is re- 
quired. In the third case, o< exists as a voltage, 
but a coefficient is needed. This means that a 
quarter- square multiplier must be used, as shown 
below. 



The circuit below meets all these requirements. 
The RESET and OPERATE times are individually 
adjustable over a wide range. The iteration is 
started and stopped by a manual switch. This 
switch allows the entire iteration to be switched 
through RESET, HOLD, and OPERATE modes. 



tt- 



I 



Figure 3—5 

The net result is that an input is multiplied by a 
constant coefficient, « /lO, which is a number be- 
tween zero and one (in magnitude). This circuit 
is, therefore, equivalent to feeding x through apot, 
except that changing the value of « will reset 
the pot" automatically. 



-20 VOLTS 



(RELAY SUPPLY) SI 



TO RESET COIL OF AOO 
TO OPERATE COIL OF AOO 



+10 



TO RESET COILS OF 
PROBLEM INTEGRATORS 



TO OPERATE COILS OF 
PROBLEM INTEGRATORS 




Finally, this circuit, which involves only one inte- 
grator, two pots, and a comparator, can be used on 
machines not equipped for repetitive operation. 



d) Mode Control . For the purposes of iterative 
computation, it is necessary to cycle the com- 
puter mode automatically between RESET and 
OPERATE. The repetitive operation unit does 
exactly this, and can be used in some forms 
of iterative computation. However, its use 
suffers from a number of limitations: 

1) For some purposes it is too fast - the 
longest OPERATE cycle obtainable is 200 
milliseconds, which makes readout on an 
X-Y plotter impractical. Also, manypro- 
grams rely heavily on relay comparators 
to make the necessary logical decisions, 
and a 200 millisecond OPERATE cycle 
followed by a 10 millisecond RESET cycle 
does not always allow these relays suffi- 
cient throwing time. 

2) When using the rep-op drive, there is 
some uncertainty about the first and last 
cycles. Switching transients may destroy 
the accuracy of the first computation. In 
**open-loop" iteration, (parameter sweep) 
this is not a serious defect, butin **closed 
loop" iteration, in which the parameter 
values for a given OPERATE cycle depend 
on the stored results from previous cycles, 
the effect can be serious. 

3) It is desirable to have a means of cycling 
the computer mode on machines not equip- 
ped with rep-op. 



When using this form of relay switching, the indi- 
vidual integrator coils should be patched directly 
to the switch and relay contacts. The regular mode 
busses are not used. 

Switch SI is initially in the LEFT position putting 
AOO in RESET and all integrators which are patched 
to the relay contacts in HOLD. 

Putting switch SI to the CENTER position puts 
every integrator in HOLD. When SI is thrown to 
the RIGHT position, AOO integrates and produces 
an asymmetrical sawtooth, as shown in Figure 
3-6, AOO starts from an initial condition + 10a 
and integrates down at the rate of -10b volts per 
second. When it reaches -10 volts, comparator 
34 switches, and AOO integrates up at the rate of 
+ lOab volts per second until it reaches + 10a, 
when comparator 34 switches again. The integrators 
patched to the comparator contacts are cycled 
between RESET and OPERATE. 



V(t) 




Figure 3-6. Output of AOO Versus Time 
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The output of AOO varies from + 10a to -10 volts. 
The total voltage variation is 10a + 10 or 10 (a+ 1) 
volts. With Comparator 34 in the position 
(which puts the computer on RESET) the rate of 
integration is -10b volts per second, and hence 
if T^ is the length of the RESET cycle the result 
is 10b = 10 (a+ 1) or 



a+ 1 



R 



Similarly, the integration rate of AOO is +10ab 
volts per second when comparator 34 is inthe 
position and the computer is in OPERATE. So, 
therefore, lOab Tq =10 (a+1) or 

a+ 1 



ab 



. , .where Tq is the length of the OPERATE cycle. 
Solving for a and b. 



T + T 
R O 

T T 
R O 



If a ten-second OPERATE cycle followed by a 
one-second RESET cycle is wanted, a =0.100 and 
b =1,100 is necessary. Since apot-setting greater 
than 1 is impossible, a gain of 10 on the integrator 
is used and the '*b" pot is set for b/10 = 0.110. 

The operational procedure is as follows: put the 
manual switch to the left, and select the HOLD 
mode on the pushbutton mode- control panel. It 
can be regarded that the L, C, and R switch posi- 
tions represent RESET, HOLD, and OPERATE 
cycles for the entire iteration ; RESET and OPER- 
ATE for the individual runs are obtained by the 
switching of comparator 34. 

In section 30 an integrator was used for automatic 
parameter adjustment. This integrator must cycle 
between HOLD and OPERATE during the iteration, 
but should be reset at the beginning of the itera- 
tion. This object can be achieved by patching its 
OPERATE coil to the + contacts and its RESET 
coil to the left contact on the manual switch. The 
integrator will then be reset with AOO. 

Since the iteration commences with a RESET cycle, 
the parameter, « , will be updated by an amount 
A «: before the first run. If is the desired value 
of « for the first run, the IC on the integrator should 
be «^ -Aoe . 
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4. ITERATION SCHEMES 



In cases where a parameter, « , is being used to 
meet a boundary value or a minimum condition, 
it must be decided what change to make in « for the 
next run. A formula for calculating the parameter 
value for the next run from the results of previous 
runs is called-an inte ration scheme . Two iteration 
schemes will be considered here. For simplicity, 
they will be limited to the case of a single param- 
eter, oc , Iteration with several parameters is con- 
sidered in References 6 and 7. 

a) Correction Proportional to Error . Suppose 
it is necessary to vary a parameter or initial 
condition, « , to meet the boundary condition 
Xf =a, where is the final value of a variable 
X, such as the horizontal displacement of the 
bullet in the missile trajectory problem or 
the temperature T in the heat equation. 

If an error is defined by f =X£ - a, then a reason- 
able change in « would be 

Aoc="G;,:: 

where G is a constant. This scheme makes large 
corrections for large errors and small corrections 
for small ones, which is reasonable. If « eventually 
settles down to a steady value, thenA«will approach 
zero, and therefore e will approach zero, which is 
what is wanted. ' V t A / \ ^ 

If e is an increasing function of «^ then the con- 
stant G should be negative so as to change « in the 
right direction to decrease e . On the other hand, 
if e is a decreasing function of « , then G should 
be positive. This means that the method is most 
useful in cases where it is known a priori whether 
€ is increasing or decreasing with « . This is the 
case with the heat equation. Figure 4-1 shows how 
Figure 2-4 may be modified by adding storage and 
parameter change circuits to effect the iteration 
scheme. The parameter oc to be varied is cx =t (0) 
and the computer effects the change in « by inte- 
grating the error e during the RESET cycle, whose 
length is T-^. 

In cases where ^ is sometimes an increasing func- 
tion of « and sometimes decreasing, the above 
scheme is less practical. Inthese cases an alterna- 
tive is to minimize the absolute error] e |. This 
approach reduces the boundary value problem to a 
problem in optimization. 

b) Fixed^Stepir Size Iteration . Consider the prob- 
lem of choosing a value for the parameter « to 
minimize the result X. X may be the absolute 




error in a boundary value problem, as describ- 
ed in Section 4a above, it may be the error in 
a curve fitting problem, as in Section 2d, or 
the amount of fuel consumed during a satellite 
maneuver, as in Section 2b. Many iteration 
schemes are possible, but the following one 
is probably the simplest to mechanize. (In 
the following paragraph, subscripts refer to 
computer runs. Thus is the value of « used 
in the first run, X^ is the result obtained from 
the first run, etc.) 

1) Pick a starting value for « , say «■£ 

2) Make a run with this value and store the 
result X-|L 

3) Change by a fixed amount A « and make 
J I ; a new run with the value =oc2 +A«» 

Store the result X2. 

4) Compare the values X^ and X2. If X2 
< X^, increase « by the same amount 
Aoc and repeat. 

5) Continue as long as X decreases with each 
successive change in « . 

6) When finally an increase in X is obtained 
this is an indication that the minimum 
value has been overshot. When this hap- 
pens, cut the step-size Aoc in half, and re- 
verse its sign. 

The following table gives an indication of what to 
expect when this procedure is followed. In this 
table, X, « , and Aoc are all computer variables 
and are given in volts. 
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Run # 


oc 


X 


A<x for next run 


1 


1.00 


+ 9,82 


+ 1.00 


2 


2.00 


+ 7.17 


+ 1.00 


3 


3,00 


+ 6.61 


+ 1.00 


4 


4.00 


+ 6.29 


+ 1.00 


5 


5.00 


+ 6.73 


-0.50 


6 


4.50 


+ 6.10 


-0.50 


7 


4.00 


+ 6.29 


+ 0.25 


8 


4.25 


+ 6.07 


+ 0.25 


9 


4.50 


+ 6.10 


-0.125 


10 


4.375 


+ 6.02 


etc. — 



Note that some of the runs are redundant; run 7 
duplicates run 4 and run 9 duplicates run 6, but 
this duplication is not a serious drawback. To 
mechanize the reversal of sign, the following cir- 
cuit can be used: 

The mode control relays on integrators 1 and 2 
are patched so that they integrate in RESET and hold 
in OPERATE . The comparator contacts assure that 
the integrators will receive inputs only when Xj> 
Xn-i; that is, only when the last parameter change 
made things worse instead of better. When X^ > 
-^n-l» integrators will integrate for one RESET 
cycle, and then HOLD. 

The system in Figure 4-2 produces a damped sine- 
wave output. The pot-settings can be chosen so 
that the half-period of this damped sinewave is 
equal to the RESET time, and the amplitude decays 
by a factor of 1/2 in this time. At the end of the 
RESET cycle, the integrators hold the new value 
of A«, which is - 1/2 times the previous valnfi. 
Hence, in a given RESET cycle, one of two things 



will happen: If ^n^^^n-l value of « will be 
updated by an amount , If^n'^'^n-l* step- 
size, A« ^ is halved in magnitude and reversed in 
direction. 




2^ Of" ^ SJ OF AMPL *i 



INTEGRATOR 

NETWORK 
FOR AMPL #1 




RESISTOR. 
NETWORK 



Figure 4-2 

Note that in this second case, the value of A« is 
changed, but the value of « is not. Hence, the next 
run will be made with the same value of « , and 
should yield the same value of X. 

The comparator will then be comparing two equal 
voltages, and it becomes a matter of chance whether 
it will decide to update « or reverse Aoc during the 
next RESET cycle. To avoid this difficulty, the 
comparator is slightly biased so that when in 
doubt, it updates «. 
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5. PROGRAM FOR ITERATIVE SOLUTION OF A BOUNDARY- VALUE PROBLEM 
USING ONLY ANALOG COMPONENTS 



In this section the trajectory problem given in 
section 2c is solved. The circuit is described in 
detail, and a typical set of results is given. 

The dynamic equations of the bullet, which were 
given in Section 2c, are repeated here for reference: 



1) The circuit was chosen to illustrate a 
number of interesting techniques. A num- 
ber of modifications could have been made 
in the interest of equipment economy, but 
the present circuit illustrates the indi- 
vidual operations better. 



X =Vo cos 9 x(o) =0 (5-1) 

y = -g y(o) =0; y(o) = Vq sin 9 (5-2) 

The solution of these equations on an analog com- 
puter is straightforward. Equation 5-1 requires one 
integrator and equation 5-2 requires two. An in- 
spection of figure 5-1 indicates that an additional 
fifteen amplifiers and three multipliers are re- 
quired to effect the iteration scheme. In view of 
the triviality of the differential equations (Eq. 5-1 
and Eq. 5-2 may be solved analytically by direct 
integration) it appears ridiculous (and indeed it is) 
to use such a sophisticated and expensive iteration 
scheme. However, the circuit is intended as a 
general illustration, not as a solution of a specific 
practical problem. The following points should be 
noted: 



2) 




The interest here is in the iteration 
scheme, not in the differential equations. 
Hence a trivial set of differential equa- 
tions, which leads to an interesting and 
complex iteration scheme, was delib- 
erately chosen. The dynamic equations 
could have been made considerably more 
complex by the addition of aerodynamic 
forces on the projectile. Some of these 
aerodynamic terms are quite complex and 
nonlinear, and would result in the bulk 
of the computing capacity being used 
to solve the differential equations. The 
iteration scheme would, of course, remain 
the same. However, from the point of view 
of illustrating the iterative techniques 
nothing would be gained by making the dif- 
ferential equations more complex. 



RESET SIGNAL FROM 
COMPARATOR 34 



•flO0C=-fiOcot^ 



505 

J:2xl0-4x 




2R/I05 



so' 




Figure 5-1. Program for iterative Solution of Boundary- Value Problem 
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3) Most of the amplifiers (and two of the three 
multipliers) are performing essentially 
digital functions — storage, gating, logic, 
etc. The use of a few digital devices (AND 
gates, flip-flops, etc.) would materially 
simplify the program. Such components 
are available for larger computers, such 
as the 100-amplifier PACE 231-R. 

The EAI DES-30 Digital Expansion System, 
a general-purpose, patchable digital logic 
package, provides digital expansion of the 
TR-48 analog computer used here. Ap- 
pendix 3 describes this essentially hybrid 
implementation of the problem in which 
the all- analog simulation is modified so 
that the logic functions are performed by 
DES-30 logic elements. 

4) In many practical cases, it is possible to 
use a far simpler iteration scheme than 
the one shown here. The iteration circuit 
in Figure 4-1, for solving the heat equa- 
tion, should be regarded as more typical 
of a large number of one-point boundary- 
value problems. In this example also, the 
differential equation would be made more 
complex if appropriate non-linearities 
were introduced— such as a temperature- 
dependent thermal conductivity. Such a 
nonlinear ity would result in the addition 
of a diode function generator and a multi- 
plier to the circuit for solving the dif- 
ferential equation, but the iteration scheme 
would remain the same. 

ANALYSIS OF THE CIRCUIT 

The parameter to be controlled is the firing angle 6. 
However, the equations require sin 6 and cos^. 
If 6 were chosen as the iteration parameter, it 
would be necessary to use two DFG's to generate 
sin 9 and cos 9. The equipment requirement is 
reduced by defining a = cos d and iterating on a. 
Then, sin d can be obtained from the relation 



"n+1 n+1 n 



sin 9 =y l-cos^g, which holds when 9 is in the 
first quadrant. This relation can be mechanized by 
a single amplifier and one multiplier (see ampli- 
fier 18, figure 5-1). This circuit makes use of the 
possibility of separating the TR-48 quarter- square 
multiplier into two independent generators. The 
iteration scheme used is proportional correction: 

"n+1 

where and €^ are the values of cc and e for the 
nth run. 



The use of the plus or minus sign is determined by 
the following rule: The signs of Aan ACn 
compared; if they are alike, this is taken as an in- 
dication that da/de is positive, and the - sign is 
used. If they are unlike, the + sign is used. 

Figure 5-1 gives the appropriate circuit. Ampli- 
fiers 06, 07, and 10 solve the dynamic equations. 
The value of x is to be sampled when y crosses 
zero. This will not happen at the end of an operate 
cycle, but at some time during the cycle. Hence a 
comparator-controlled track- store imit is used. 
When y crosses zero, comparator 35 flips to the 
minus side, switching amp 03 from TRACK to 
STORE, Amplifier 07 continues to integrate down- 
ward (this means a negative value of y, which has 
no physical significance). A feedback diode is used 
to prevent the output of amp 07 from becoming 
more negative than about half a volt. This is to 
prevent overload in case the OPERATE cycle con- 
tinues for a long time beyond the collision point 
(the length of the OPERATE cycle is fixed) . 

Amplifier 03 is a combination summer and track- 
store unit (see Appendix 4) . The desired value of 
the range R is set in on pot 00, and amplifier 03 
tracks and stores the error. At some point in the 
OPERATE cycle amplifier 03 is switched from 
the TRACK mode to the STORE mode, and it 
should remain in this mode throughout the rest of 
the OPERATE cycle and the following RESET 
cycle. However, during this RESET cycle, the con- 
trol signal, which is derived from amplifier 07, 
will be reset to zero, which will flip the com- 
parator to and return amplifier 03 to the 
TRACK mode prematurely. To prevent this, the 
second input to comparator 35 is connected to the 
RESET signal from comparator 34 (see figure 5-2) 
which is performing the mode cycling. When the 
main problem integrators 06, 07, and 10 are being 
reset, a -20 volt signal appears at this terminal, 
latching the comparator in the state and as- 
suring that amplifier 03 remains in the STORE 
mode throughout the RESET cycle. Since the -20 
volt relay power supply is a larger signal than 
is normally used on the TR-48 patch panel, an 
input resistor of lOK is used to cut down com- 
parator sensitivity. This resistor may be a **gain 
10" resistor from any of the unused amplifier 
networks. 
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Amplifiers 14 and 22 delay the stored error by one 
cycle (see Table 3-1 and Figure 3-2). Amplifier 14 
(marked S*) behaves like an ordinary S (Storage) 
amplifier, but its mode control is a bit different 
from that described in section 3a. It cycles be- 
tween OPERATE and HOLD instead of OPERATE 
and RESET. An ordinary Storage amplifier, as de- 
scribed in section 3a, could have been used, but 
this trick enables the amplifier to share a dual 
integrator network with amplifier 15, reducing the 
number of integrator networks required. 

Amplifier 08 is a high-gain amplifier. Instead of a 
feedback resistor, it has a feedback limiter circuit 
designed to limit it to approximately ±10 volts. The 
output of amp 08 may be thought of as essentially a 
digital signal, having the value -10 when € ^-^ ^n-1 
and +10 when < 1, or in more concise 
terms, its output is -10 sign (ACn). Similarly, the 
output of amplifier 16 is +10 sign Aa^j, which has 
been computed from the value of La for the last 
run, delayed one cycle. The signs of Aq:^ ^^n 
are compared by multiplier 13, which acts as a 
''coincidence gate" or ''exclusive OR" circuit. 
The output of amplifier 13 is -10 volts if Aa^ and 
ACn have the same sign, and +10 volts if they 
have opposite signs. Multiplier 21 multiplies this 
logic signal by c^, so that the output of amplifier 
21 is proportional to ±en, the sign depending on 
the relative signs oi^OL^ and ACn. The signs on 
the multipliers are chosen so that the output is 
proportional to -Ao^n+l' '^^^^ signal is delayed 
by amplifiers 11 and 23 and converted to a binary 
signal by amplifier 16, so that it can be used to 
determine the direction of change for the next run. 



The output of amplifier 21 is just what is needed 
for the iteration scheme i^otxi+l ~ ~ "^^n 

SIGN (Ao^jj ^^n'- introduces the constant 

G and the appropriate scale factor. Its output is 
-10 A^n+l/TR. Amplifier 15, the parameter up- 
dating integrator, is connected so that it operates 
when the main problem is in RESET and holds 
when the main problem is in the OPERATE mode. 
It integrates the constant -10 A(%+i/Tj^ for a con- 
stant interval Tj^ (the RESET time) and its output 
changes by an amount +10 Ao^n+i in this time. Thus 
its output, +10 a = +10 cos 9 is appropriately up- 
dated during the RESET cycle and its new value is 
held during the following OPERATE cycle. Multi- 
plier 25 and amplifier 18, as described above, pro- 
vide the necessary squaring, subtraction, and 
square root functions to generate 10 sinB from the 
input 10 cos ^. 

Figure 5-2 gives the circuit for the mode control 
cycling. The operation of this circuit was analyzed 
in detail in section 3-d. Note that the MASTER 



Tr/To 



+ 10 




34 L 



SWITCH JR 
-20 NO. I 



RELAY 
VOLTS 



TO OPERATE COILS OF MAIN 
►PROBLEM INTEGRATORS 

>T0 RESET COILS OF MAIN 
PROBLEM INTEGRATORS 



►TO OPERATE COIL OF A02 

^TO RESET COIL 0FA02. THIS IS 
MASTER RESET TERMINAL. CONNECT 
TO INTEGRATORS THAT MUST BE 
RESET WITH THE ENTIRE INTEGRATION 
SCHEME. 



Figure 5-2. Mode Control Circuitry for Boundary-Value 
Problem 



RESET terminal on switch #1 may be patched to 
the RESET coil of the parameter-updating inte- 
grator 15, to enable a first-run value a-^ to be 
established on the integrator. The computer mode 
busses are completely unused, and all mode con- 
trol signals are obtained from the manual switch 
and the comparator. The -20 volt relay power 
supply is terminated on the TR-48 patch panel in a 
terminal marked RV (relay voltage). This terminal 
is in the lower left corner of the readout panel. 

Since the RESET and OPERATE busses are not used 
in this circuit, the operation of the circuit will be 
independent of the mode pushbuttons on the control 
panel, except for the POT-SET mode, which closes 
a relay around every amplifier. The POT- SET mode 
is therefore a convenient one for setting all storage 
amplifiers to zero before the start of the iteration. 
The sequence of operations to be performed by the 
operator is as follows: 

1) Put the computer in the POT-SET mode to 
clear all storage to zero. 

2) Put SWITCH #1 to left position. 

3) Depress the HOLD pushbutton to bringthe 
computer out of the POT-SET mode. 

4) Put SWITCH #1 to the right position. 

The computer will now cycle back and forth between 
RESET and OPERATE, updating the value of the 
parameter a each time. The results may be re- 
corded on a X-Y plotter, and the convergence 
observed. 
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Results for a typical series of runs are given in 
figure 5-3. The target was placed at 2,000 feet, 
and a first-run value of a was chosen that gives 
an extremely large initial overshoot. Note that the 
results of successive runs approach the target mon- 
otonically, and that the amount of the correction 
decreases as the larget is approached. 

Table 5-1. Potentiometer Settings 



Pot 
Number 



Parameter 
Description 



Numerical 
Setting 



It is worthwhile to pay some attention to the mean- 
ing of the constant G. The entire iteration scheme 
may be thought of as a feedback control system con- 
trolling the parameter o^. The correction is pro- 
portional to the error e , and the constant G may 
be thought of as the gain of the iteration scheme. 
It is interesting to run the iteration scheme several 
times and observe the effect of the gain G on the 
convergence. As expected, low gain leads to slow 
convergence. Increasing the gain speeds the con- 
vergence, up to a point, but excessive gain produces 
divergence. 



00 
01 
02 
03 
04 
05 
06 
07 
08 
09 
10 
14 
19 
24 
29 



2R/10 

g/ioV^ 

1/50/5 



Limiter 



2V 



J- 



T^ /T 
R/ op 

^ R op' 
Limiter 

50,000G/T 

Limiter 

Limiter 

Limiter 

Limiter 

9,000 



r/ R op 



R 



0.040 
0.161 
0.100 
0.750 
0.470 
0.100 
0.075 
0.500 
0.300 
0.470 
0.040 
0.470 
0.470 
0.470 
0.470 



Table 5-2. Summary of Integrator Mode Patching 
02, 03 From SWITCH 1. 



06, 
10, 



07 
11 



14. 15 



21. 23 



From lower contacts on comparator 34 
RESET coils from + terminal 
OPERATE coils from - terminal 

RESET coil patched to SWITCH 1, left 
terminal. 

OPERATE coil from + terminal, com- 
parator 35 

From lower contacts on comparator 34 

RESET coils from - terminal. 

OPERATE coils from + terminal. 
(The mode of these integrators is com- 
plementary to the mode of the main 
problem integrators) 



COMPUTER RESULTS FOR 6 ITERATIONS OF 
BOUNDARY-VALUE PROBLEM, USING PRO- 
PORTIONAL CORRECTION. TARGET PLACED 
2,000 FEET FROM LAUNCHING SITE. 
NUMBERS ON CURVES REFER TO SUCCESSIVE 
.OPERATE CYCLES 




0 1,000 2,000 3,000 4,000 5,000 6,000 7,000 8,000 9,000 10,000 11,000 12,000 13,000 14,000 

X (FEET) 



Figure 5-3. Results for a Typical Series of Runs 
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APPENDIX 1 



PARAMETER SWEEP CIRCUIT 



The NORMAL/SPECIAL option on the TR-48 inte- 
grator enables us to construct a very simple cir- 
cuit for the parameter sweepproblem in Section 3 a. 
In Figure A2-1, integrator #1 is a SPECIAL inte- 
grator, while all the others, including the T and S 
(track and store) integrators are operating in 
rep-op. 



TO ARM OF 
X-Y PLOTTER 



50 



SYSTEM OF 
DIFFERENTIAL 
EQUATIONS 
(RER-OR) 



X(t) 



I TO PEN 
SQ>^0F X-Y 



PLOTTER 



Thus, the rep-op circuit will solve the differential 
equations of the system for a virtually constant 
value of oc , and the S-amplifier will store the final 
value X. After the next solution, the slightly dif- 
ferent value of will cause a small change in the 
value of X, and the output of the S-amplifier will 
change accordingly. With the « r- integrator con- 
nected to the arm and the S-amplifier connected to 
the pen, the X-Y plotter will produce a graph of 
the final value X versus the parameter « . 

The graph should look something like Figure A2-2; 
that is, the graph will not be smooth, but will con- 
sist of many small steps. However, if « varies 
over the desired range in 50 seconds, and the rep-op 
rate is 20 solutions per second, the graph should 
consist of 1,000 steps, and they should be so small 
as to be barely noticeable. 



Figure A2-1 

After every rep-op cycle, the S-amplifier stores 
the final value of the variable X. This final value 
will change gradually as changes. Since « is 
changing very slowly (it takes 50 seconds to go 
from zero to + 10 volts) it will appear very nearly 
constant from the point of view of the rep-op cir- 
cuit. In particular, if the rep-op rate is 20 solu- 
tions per second, « will change by only 0.1 volt 
during a single computation. This is the smallest 
change that can be detected on the TR-48. 



X FINAL 



Figure A2-2 
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APPENDIX 2 



TR-48 RELAY MODE CONTROL 



The mode control on the RMC TR-48 is accom- 
plished with relays which are energized from a 
-20 volt DC power supply. The individual relay 
coils appear on the patch panel, which allows for 
independent control of integrator mode. 

The relays normally receive inputs from two mode 
control busses, called the RESET bus and OP- 
ERATE bus. Each bus is energized from the corre- 
sponding pushbutton on the mode control panel: 

1) When the RESET mode is selected with a 
pushbutton, the RESET bus is energized 
with -20 volts and the OPERATE bus is 
de-energized. 

2) When the OPERATE mode is selected, the 
OPERATE bus is energized, and the RE- 
set bus de-energized. 

3) In the REP-OP mode, the RESET and 
OPERATE busses are energized alter- 
nately, thus cycling the computer between 
RESET and OPERATE. 

4) When the HOLD mode is selected, both 
busses are de-energized. 

Each dual integrator network has two input termi- 
nals on the patch panel which go to relay coils. 
These terminals will be referred to as the RESET 
coil and the OPERATE coil. (This is an over 
simplification — actually, one of the terminals 
energizes two relay coils, but this is not important 
from a programmer's point of view): 

1) When the RESET coil is energized, the 
integrator is in the RESET mode. 

2) When the OPERATE coil is energized, the 
integrator is in the OPERATE mode. 

3) When neither coil is energized, the inte- 
grator is in HOLD. 

4) If both coils are energized simultaneously, 
the integrator may be in any of several 
modes. Because of this uncertainty, the 
programmer should avoid simultaneously 
energizing both coils. 

In normal operation, each coil is patched to the 
corresponding bus. The mode of the integrator 



is then selected by the pushbuttons in the normal 
manner. This standard mode control connection 
can be made with a double horizontal bottle plug. 
If this plug is removed, the coils can be patched 
in many different ways. For example, reversing 
the connections (RESET coil to OPERATE bus and 
OPERATE coil to RESET bus) gives a tracking unit, 
as described in Section 3a. Integrator coils may 
also be energized directly from the -20 volt relay 
power supply or through comparators or manual 
switches. The relay power supply appears on the 
patch panel in the lower left-hand terminal of the 
readout panel, at the terminal marked RV (relay 
volts). 

For some applications, it is desirable to have 
some integrators running in Rep-op and some in 
real time simultaneously. This result can be 
achieved on the TR-48 as follows: a shaded por- 
tion of the integrator network contains two holes 
marked NORMAL and two holes marked SPECIAL. 

When the holes marked SPECIAL are patched 
together, the integrator operates normally. (When 
a bottle plug is used for this connection, the word 
SPECIAL is covered up. Hence the operator sees 
the word NORMAL when the integrator is operat- 
ing normally). When the holes marked NORMAL 
are covered with a bottle plug and the word 
SPECIAL is exposed, the integrator operates in a 
special manner, i.e., it goes into the OPERATE 
mode as soon as the REP-OP mode button is de- 
pressed, and remains in OPERATE (with a real- 
time capacitor, not a rep-op capacitor) as long as 
the computer is in repetitive operation. An applica- 
tion of such an integrator is given in Appendix 1. 

The mode control relays on the TR-48 are shared, 
by the two integrators in a dual integrator network. 
Hence the two integrators in a dual network must 
be in the same mode at any given time. The pro- 
grammer must keep this in mind when assigning 
numbers to various components. Occasionally an 
integrator network will be wasted (if an odd number 
of integrators with a special nonstandard mode 
logic are required, one extra network will have 
this same logic, whether it is needed or not). How- 
ever, the waste will never amount to more than a 
few integrator networks. Note that an S- amplifier 
(storage amplifier) has the same mode logic as a 
conventional integrator. Hence if three S- amplifiers 
are needed in an iteration, the odd network can be 
used in the problem itself. 
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One other feature of the TR-48 integrator network 
has been used in this article. This is the possibility 
of changing the feedback capacitor on an integrator 
to increase its integration rate. This feature allows 
a wider range of time- scale factors than would 
otherwise be the case. A bottle plug, called the 
j8 plug, is normally patched into each integrator. 
Removing this plug replaces the feedback capacitor 
by a capacitor of 1/10 the value, increasing the 
integration rate by a factor of 10. This feature 
operates independently of mode control patching, 
rep-op, NORMAL/SPECIAL connections or any- 
thing else. Whatever feedback capacitor the inte- 
grator would have with the jS plug (a rep-op or 
real-time' ' capacitor) is replaced by a capacitor 
of 1/10 the value. In contrast to the mode control 
connections, the dual integrator has two separate 
P plugs — one for each integrator in the network. 



Note: The TR-48 with Electronic Mode Control (EMC) differs only 
slightly from the RMC Computer in so far as the programmer is con- 
cerned, despite its completely different system of controlling the 
modes of integrators. The electronic switches, which replace the 
mode'control relays, are operated front logic signals, i.e., on the 
TR-48, a logical 7 = +5V and a logical 0 = OV. Hence, the busses 
which feed the mode switches operate at these two levels just as 
the RMC busses operate at -20V and OV levels. 



The following table illustrates EMC operation. 



INPUTS 


RESULTING 


IC 


OP 


MODE 


1 


0 


RESET 


0 


7 


OPERATE 


0 


0 


HOLD 


7 


7 


NOT ALLOWED 



These inputs can be obtained by depressing the proper mode control 
buttons if the mode-control bottle plugs are patched, or by patching 
from devices which produce these levels (electronic comparators, 
gates, flip-flops, etc.). There is a terminal in the readout panel 
labeled +5V which can be used as a logical 7 level "supply**. 

All other components on TR'48 dual integrator networKs apply to the 
EMC network as well as the RMC network. 

In addition, EMC TR'48*s are usually available with Electronic 
Comparator modules, which comprise one electronic comparator net' 
work, two digitally-controlled analog switch networks, and two track' 
store amplifier networks. The latter obviate the need for using 
standard integrators as track'Store devices and thus add to the 
flexibility of the computer. Electronic Comparator modules are also 
available for the RMC computer. 

Symbols 

7. Electronic Comparator (EC) 



ANALOG 
INPUTS 



62 



I DIGITAL 
OUTPUTS 



When ej + 62 > Q, the "V output is at the "7" level 
When ej + 62 < 0, the "7" output is at the "0" level 



2. D/A Switch 



ANALOG 
INPUT 

DIGITAL- 
COMMAND 




ANALOG 
OUTPUT 



3. T/S Amplifier 



ANALOG 
INPUT' 











T 




SAMPLED 
ANALOG 
OUTPUT 



DIGITAL COMMAND 
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APPENDIX 3 



TR-48/DES-30 HYBRID IMPLEMENTATION 
OF 

ONE PARAMETER OPTIMIZER (BOUNDARY VALUE) PROBLEM 



The all-analog simulation of the trajectory prob- 
lem given in Section 2c is modified for hybrid 
implementation. The necessary logic functions are 
performed by the logic elements of the EAI DES 
-30 Digital Expansion System 

ANALOG PORTION 

The analog portion of this problem is the part of 
figure 5-1 used to solve the equations 5-1 and 5-2. 
It is repeated here for reference. 

MODE CONTROL 

The mode control program (see Figure A3- 2) 
consists of DC A and FF C. Its purpose is to 
provide signals to alternately throw the TR-48 into 



OP (9 seconds) and RS (9 seconds). The 9 sec 
period is dictated by the problem solution time — 
longer periods may be used if desired, but not 
shorter periods. The program also provides the 
signal necessary to control the OP-HD mode of 
integrator 15 — the one that operates a = cos 9 
from AQin+ 1 ^i^^ii^g the time that the problem is in 
RS, and goes to HD when the problem is in OP. 
See the timing diagram in Figure 3A-3. 

The LP permits manual start- stop. When the LP 
is pushed, DC A counts down at the rate of 1 count 
per second, producing a Cq blip and reloading 
itself to 9. The Cq blip triggers FF C so that its 
output is high for 9 seconds and low for another 9 
seconds. Thus, the complementary outputs of FF C 
are just the signals necessary to control the OP, 
RS mode inputs. 




Aan+I, FROM D/A AMR B 



OM 
DES-30 



FROM 
DES-30 



Figure A3-1: Analog Program for Equations 5-1 and 5-2 
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OPTIMIZATIQN LOGIC 

T/S A tracks the error, e^, during the OP mode, 
while the altitude y is still positive (above ground) . 
When y goes negative (the projectile hits the ground) , 
comparator A goes low, and T/S unit A holds. 



OP-RS 
TIMING 



OPp 
RS I 



-€n TRACK p 

(T/S A) STORE I 



•9 SEC- 



r 



After the OP period is finished, the RS mode 
output throws T/S B into TRACK. It then tracks 
the output of T/S A (with sign inversion) . Thus, on 
the next OP cycle, T/S A tracks the present error, 
Cn, while T/S B holds the old value, ^n-l- See 
timing diagram in Figure 3 A- 3. 

Comparator B goes high whenever ^^-1 " ^n ^ 
(that is, when A c is positive) and remains low when 
Ae is negative. Thus, its output may properly be 
called **sgn Ae". Flip flop B stores this informa- 
tion as of the instant when T/S A goes to STORE — 
that is, at the earliest instant both T/S A and T/S 
B are storing. This is achieved by disabling this 



fn;l ^ TRACK 
(T/S B) 



STORE 



HI 



SAMPLE 
SGN 



SAMPLE 
SGN AQn 



INTEG 
15 



_r 



HO p 
OP I 



l_ 

r 



r 



Figure 3A— 3: Timing 
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FF when T/S A goes to STORE. The same effect 
can be achieved by enabling FF B by a blip ob- 
tained from differentiating the trailing edge of 
Comp. A output, but this would require a differen- 
tiator. 

Note that FF B * 'tracks" sgn Ac throughout the 
time when y> 0, It might be asked whether this 
has any adverse effect. It can be seen that this is 
not the case because integrator 15 is in HD through- 
out the period y>0, so that it doesn't care about 
sgn Ac (or anything else, for that matter). 

FF A is e nabled by a blip which is the AND of Cq 
and FF C — that is, a blip occurring just prior to 
going into the OP mode. At this time, the output of 
amplifier B has Ao^n on it, comparator C output is 
sgn Ao^, and so FF A stores sgn Aa^, Recall that 
Aq^ is the quantity used to generate the parameter 
of interest, o^i^^i, that is, Ao^^ is determined 
during the OP mode of the previous run, not the one 
we are just going into. 

The outputs of flip flops A and B are compared in 
a negated XOR (Exclusive-Or) program. An XOR 
is basically a digital comparator — its output is 
high whenever the two inputs are unequal (one is 
low while the other is high). Thus the negated 
XOR output is high whenever the inputs are equal — 
both high or both low. That is, the signal 6 in 
Figure 3A-2 is high whenever sgn Ae^and sgn 
are equal; in Boolean notation 6 = (sgn Ae^^). 
(sgn AcXj^) + (sgnA^j^) (sgnAa^)? algebraic 
notation 6 = (sgnAe^) • (sgnAa^)- Thus, the XOR 
programs replaces multiplier 13 in the all-analog 
implementation of this program. 

Amplifier B and the associated electronic switches 
control the selection of the sign in the equation 
Ao^+ 1 = ± Ge in accordance with the signal 6 
just described: when 6 =1, indicating that sgn A^n 
and sgn Aa^^ are equai, the output of amplifier B 
is when 6=0, -6^. (Recall that integrator 15 
requires -Ao^^^^ at its input). Thus, amplifier B 
and its switches replaces multiplier 21 in the all- 
analog simulation. 



INITIALIZATION 

In all programs of this type, initialization should 
carefully be considered. Initialization consists of 
two problems: 1) How does one get the problem 
started? (This is no more a simple analog set-up; 
we also have a digital section to worry about, 
and simply pushing the OP button on the TR-48 



won't do!) 2) How does one make sure that all 
initial conditions are proper, both analog and 
digital? 

The question of starting (and, incidentally, the 
related question of stopping) has been taken care 
of by the LP in the mode control section. When the 
LP is released, FF C is cleared, and the TR-48 is 
in RS, This assumes that the clock is running and 
the TR-48 is in SLAVE — which is the normal 
state of affairs in hybrid simulations when running. 
When the LP is latched, DDC A starts counting and 
the simulation is under way. 

The question of digital initial conditions is that of 
the correct initial states of flip flops (GPR FF, 
DC FF, MT & DIF FF, etc.) For example, we may 
wish to set DC A to 09 manually so the cycle 
starts in the RS mode. If DC A is cleared, the 
cycle starts with OP mode. But most importantly, 
we must decide whether FF's A and B should be 
cleared or set to begin with. If we know what sgn 
A« ^ and sgn Ae ^ are for the first two run (and in 
this problem we can find this out because it hap- 
pens to be a very simple one), then we can set 
FF's A and B accordingly. If we do not know, as 
will be the case in a more sophisticated problem, 
then we must program an initialization run to 
determine the initial quantities. In this particularly 
simple situation, we can ignore the initial setting of 
FF's A and B altogether. This results in at most 
2 bad runs (the first two) in which the optimization 
may go in the **wrong direction". However, after 
T/S A and T/S B have had a chance to operate 
once, things will quiet down. 

A LIST OF MNEMONICS (DES-30) 



c 

o 


Carry Out 


DC 


Down Counter 


DIF 


Differentiator 


EMC 


Electronic Mode Control 


FF 


Flip Flop 


GPR 


General Purpose Register 


HD 


Hold Mode 


LP 


Latching Pushbutton 


MP 


Momentary Pushbutton 


MT 


Monostable Timer 


OP 


Operate Mode 


RS 


Reset Mode (IC) 


T/S 


Track/Store 
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APPENDIX 4: THE COMBINATION SUMMER 
AND TRACK-STORE UNIT 



In many applications, it is necessary to store the 
sum of two or more signals. This can, of course, 
be accomplished with two amplifiers - a summer, 
and a track-store unit. However, the functions of 
both may be combined in a single track- store unit 
with several input resistors (Figure A4-1). 




Figure A4-1 

When the switch is closed, we have 



y = - 



1+RCS 



Ar^ R. i 



. , .and the unit acts as a combination track- 
summer. 

The point marked SJ"^ in Figure A4-1 appears on the 
patch panel area of a TR-48 integrator network so 



that additional IC input resistors can be patched in 
to take advantage of this capability. In many cases, 
this feature saves one or more amplifiers. 

Since the input resistor network of a track- store 
integrator is not normally used, it can be patched 
to the SJ'^ and used for multiple inputs. A symbol 
for this type of patching is given in Figure A4-2. 



NUMBER OF SUMMING 
^NETWORK USED. 

-^14 




Figure A4-2 



Note that a '*gain of 10" on a TR-48 means a lOK 
input resistor, which produces an actual gain of one 
in this case, since the feedback resistor is also 
lOK. (Some TR-48' s have 50K input and feedback 
resistors in the integrator IC networks. The in- 
dividual computer should be checked for this before 
the circuit is used.) 
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